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ABSTRACT 

Context. The role played by protostellar feedback in clustered star formation is still a matter of debate. In particular, protostellar 
outflows have been proposed as a source of turbulence in cluster-forming clumps, which may provide support against global collapse 
for several free-fall times. 

Aims. Here, we seek to test the above hypothesis in the case of the well-documented NGC 2264-C protocluster, by quantifying the 
amount of turbulence and support injected in the surrounding medium by protostellar outflows. 

Methods. Using the HERA heterodyne array on the IRAM 30m telescope, we carried out an extensive mapping of NGC 2264-C in 
the three molecular line transitions '2CO(2~l), "C0(2-l), and C"*0(2-l). 

Results. We found widespread high-velocity '^CO emission, testifying to the presence of eleven outflow lobes, closely linked to the 
compact millimeter continuum sources previously detected in the protocluster. We carried out a detailed analysis of the dynamical 
parameters of these outflows, including a quantitative evaluation of the overall momentum flux injected in the cluster-forming clump. 
These dynamical parameters were compared to the gravitational and turbulent properties of the clump. 

Conclusions. We show that the population of protostellar outflows identified in NGC 2264-C are likely to contribute a significant 
fraction of the observed turbulence but cannot efficiently support the protocluster against global collapse. Gravity appears to largely 
dominate the dynamics of the NGC 2264-C clump at the present time. It is however possible that an increase in the star formation rate 
during the further evolution of the protocluster will trigger sufficient outflows to finally halt the contraction of the cloud. 

Key words, stars : formation : circumstellar matter - ISM : individual objects : NGC 2264-C - ISM : kinematics and dynamics - 
ISM : outflows 



1. Introduction 

1.1. Clustered star formation, turbulence and protostellar 
feedback 

It is now well established that a large fraction of young stars 
in giant molecular clouds form in groups and clusters rather than 
in is olation (e.g. |Carpenter||2000[ [Adams & Myers||2001[ |Lada 
[&Lada 2003i ). 

Three main classes of models have been proposed for reg- 
ulating the star formation process and explaining the origin of 
the stellar initial mass function (IMF) in young star clusters. The 
first scenaiio is based on turbulent fragmentation of the parent 
molecular cloud (e.g. ,Padoan & Nordlund|2002 Hennebelle & 



Chabrier||2008 1, which produces an IMF-Uke core mass distri- 
bution as observed in the nearest cluster-forming regions (e.g. 
Motte, Andre & Neri| [T998] [S tanke et al. 2006). Briefly, self- 
gravitating pre-stellar condensations (each containing one lo- 
cal Jeans mass) form as turbulence-generated density fluctua- 

* Based on observations carried out with the IRAM 30m telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany), and 
IGN (Spain). 

** '^C0(2-l) data cubes in FITS format are available at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
[http://cdsweb.u-strasbg.fr^ 



tions, then decouple from their turbulent environment, and even- 
tually collapse with little interaction with their surroundings. In 
this picture, the IMF results primarily from the properties (e.g. 
power spectrum) of interstellar turbulence, which also controls 
the global star formation rate and efficiency (e.g. Mac Low & 
Klessen 2004, |Krumholz & Tan|2007l ). The second class of mod- 
els emphasizes the role of protostellar turbulence and feedback 
in regulating the star formation process (e.g. Norman & Silk 



|T980t [Adams & Fatuzzo 1996 S huet al.|20 04'). Here, the IMF 
and the star formation efficiency are determined by the stars 
themselves through the collective effects of their feedback on 
both individual cores and the parent cloud. A third scenario ex- 
ists, however, according to which turbulence is only responsible 
for forming protostellar seeds but plays no direct role in regu- 
lating the star formation process and shaping the IMF, which is 
mainly determined by competitive accretion and dynamical in- 
teractions between individual cluster members (e.g. Bonnell et 
al. 2001 ; Bate et al. 2003 ). 

From an observational perspective, millimeter-wave obser- 
vations of molecular clouds have revealed supersonic linewidths, 
which are presumably due to turbulent motions. Theory suggests 
that turbulent motions can be treated as an additional pressure 
(e.g. Chandrasekhar 1951 Von Weizsacker 1951| l, so that su- 
personic turbulence increases the effective Jeans mass supported 
against collapse. 
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Protostellar outflows have long been considered a plausible 
way of driving a significant fraction of molecular cloud turbu- 
lence, especially in cluster-forming clumps, where stars form 
more efficiently than in the bulk of molecular gas ([Bally et al.| 



T996l|Reipurth & BaUy|200T)|Knee & Sandell|2000| i. Recently 
Nakamura & Li ( 2007[ l discussed the possible effects of pro- 
tostellar outflows on cluster formation. In particular, they ar- 
gued that, due to its short decay time (e.g. Mac Low et al. 



1998 1, the "interstellar turbulence" initially present in a cluster- 
forming clump is quickly replaced by turbulent motions gener- 
ated by protostellar outflows. According to [Nakamura & Li ] the 
protostellar outflow-driven turbulence dominates for most of a 
protocluster's lifetime and acts to maintain the cluster-forming 
region close to overall virial equilibrium for several dynami- 
cal times, avoiding global free-fall collapse. The exact mecha- 
nism by which outflows inject turbulence in the parent cloud is 
not well understood yet and may involved fossil cavities (e.g. 
Quillen et al.|2005 Cunningham et al.|2006j . On the scale of an 
entire giant molecular cloud, the shocks produced by outflows 
from young stars may not inject momentum and energy at a high 
enough rate to counter the rate at which turbulent energy decays 
(see Banerjee et al.[[2007] l. On the other hand, intense outflow 
activity within a young protocluster may be sufficient to support 
the parsec-scale parent gas clump, thereby modifying the out- 
come of the star formation process. 

As the role of protostellar feedback in cluster-forming clouds 
is still a matter of debate, detailed studies of the dynamical ef- 
fects of protostellar outflows in young protoclusters, where out- 
flows are particularly strong and numerous, are required to fully 
understand the process of clustered star formation. In this paper, 
we test some of the above ideas by quantifying the impact of pro- 
tostellar outflows on the dynamics of the NGC 2264-C clump, a 
protocluster at a very early stage of evolution. 

1.2. Our target region : NGC2264-C 

The NGC 2264 cluster-forming region is located in the 
Mon OBI giant molecular cloud complex at a distance d ~ 
800 pc. The cloud associated with NGC 2264 has been the tar- 
get of several millimeter line studies, including an unbiased CO 
(J=1^0) survey for molecular outflows ( jMargulis et al. 1988| l 



and a search for dense gas via multitransitional CS observations 
( [Wolf-Chase eraL[[T9^ [Wolf-Chase & Walker[[T995l l. These 
early studies revealed two prominent molecular clumps, named 
NGC 2264-C and NGC 2264-D, each associated with a CO out- 
flow. The NGC 2264 region also harbors a young star cluster 
of more than 360 near-IR sources ( [Lada et aL][1993[ [Lada &[ 
Lada|200"3] l. The most luminous object of the cluster, hereafter 



caUed IRSl after Aflen (1972), is embedded in NGC 2264-C 
and associated with the Class I IRAS source IRAS 06384-H0932 
(Lbai ~ 2300 Lq - Margulis et al. 1989). 

The internal structure of NGC 2264-C and NGC 2264-D 
was subsequently resolved by higher-resolution milUmeter and 
submillimeter studies (| Ward-Thompson et al.|2000 WilUams & 



[Garland[[2002l [Peretto et al.|2006[|2007|l. In particula r, flirough 
1.2 mm continuum observations, [Peretto et ar| (!2006) identified 



12 compact dense cores in NGC 2264-C (see Fig{T]), of which 
at least 8 are Class 0-like objects with associated near-IR H2 
jets ( Wang et al.|2002| i or compact CS outflows ( [Schreyer et al. 
'2003| l. These cores have typical diameters ~ 0.04 pc, masses 
ranging from ~ 2 to ~ 40 Mq, and column densities ranging 
from ~ 4 X 10^^ to ~ 6 x 10^^ cm 2. Based on HCO+, CS, and 



NaH^ mapping observations, [Peretto et al. ( 2006 1 also showed 
the presence of large-scale collapse motions converging onto the 



most massive core (C-MM3 with M ~ 40 M©), near the center 
of NGC 2264 -C. The IRAM Plateau de Bure interferometer ob- 
servations of |Peretto et al.[ ( [2007 ) revealed a thirteenth source, 
C-MM13, located only 13.5" (or ~ 0.05 pc) away from C-MM3. 
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Fig. 1. (a) Dust continuum map of NGC 2264-C taken at 1 .2 mm 
by [Peretto et aL] p.006} with the MAMBO bolometer array on 
the IRAM 30m telescope. Twelve millimeter peaks were found 
and interpreted as candidate Class objects (see Peretto et al. 
for details), (b) Total-power 450 fim dust continuum map of 
the NGC 2264-C protocluster taken by us with the P-ArTeMiS 
bolometer camera' on the APEX telescope. The HPBW reso- 
lution is ~ 10". The contour levels go from 1.75 to 7 Jy/10"- 
beam by 1.75 Jy/10"-beam and from 10 to 22 Jy/10"-beam by 3 
Jy/10"-beam. 



In their study, Peretto et al. (2007 1 sucessfully modeled 



the dynamical evolution of the NGC 2264-C protocluster. They 



' P-ArTeMiS is a 16 x 16-pixel prototype of the large-format submil- 
limeter bolometer camera currently being built by CEA for the APEX 
telescope (cf. Andre et al. ( 2008 1; Talvard et al. (2008 )). 



http : //www . apex- telescope . org/ instruments/pi/artemis/ 
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compared the available millimeter observations with SPH nu- 
merical simulations of the collapse of a Jeans-unstable, prolate 
dense clump. Peretto et al. ( ,200 7 ) obtained a good quantitative 
agreement with observations when starting from a low level of 
initial turbulent energy. However, the total mass of dense gas 
(i.e., with > 10"* cm"^) was ~ 10 times lower in the best-fit 
simulations than in the actual NGC 2264-C clump. When Peretto] 
et al. (2007 1 increased the mass of dense gas by a factor of 10 in 
their numerical simulations, a larger number of cores was gen- 
erated, leading to a poor match to the observations. It seemed 
therefore that an additional source of support against gravity, not 
included in the simulations of Peretto et al.'('2007), had to be in- 
voked to explain the current dynamical state NGC 2264-C. 
Peretto et al. suggested that this extra support could arise from 
protostellar feedback and/or magnetic fields. 

Because it is well documented and known to exhibit outflow 



activity ( [Margulis etar.||1988j , the NGC 2264-C protocluster is 
an ideal laboratory for probing the initial conditions of clustered 
star formation and evaluating the impact of outflow feedback 
on early protocluster evolution. We thus initiated a mapping 
study of the outflow already detected by Mar gulis et al.| ( T988| l 
in NGC 2264-C, with higher angular resolution and better 
sensitivity. Our goal was to assess the momentum injection 
rate due to outflows in this protocluster and examine whether 
outflows could affect the global dynamical evolution of the 
protocluster. 



strong source NGC7538. The observations were carried out in 
single sideband mode with an image rejection of ~ 10 dB, pro- 
viding a calibration accuracy better than ~ 10%. At 230 GHz, 
the main beam efficiency and the forward efliciency of the tele- 
scope are Befi = 0.90 and Fefi = 0.52 respectively. Throughout 
this paper, the line intensity scale adopted is in units of T^, the 
equivalent antenna temperature above the atmosphere. 



Table 1. Parameters of the IRAM 30m observations 



Line 




Frequency 


Resolution 


Bandwidth 


Mean r.m.s. 












noise" 






(GHz) 


(km.s-') 


(km.s-') 


(K) 


'2CO(2- 


-1) 


230.5379 


0.84 


157 


0.1 


'3CO(2~ 


-1) 


220.3986 


0.10 


78 


0.2 


C"*0(2- 


-1) 


219.5603 


0.10 


78 


0.3 



rms noise per spectral resolution element. 



We also mapped the NGC 2264-C region in the 450 fim dust 
continuum at ~ 10" resolution (HPBW) with the P-ArTeMiS 
bolometer array on the Atacama Pathfinder Experiment (APEX) 
telescope in November 2007. The resulting map is shown in 
Fig. [ijb)' details about the P-ArTeMiS camera and the 
November 2007 run at APEX can be found in I Andre et al.l 
( |2008| . 



2. Observations and data reduction 

2.1. Observations 

Observations of the i2cO(2-l), C"*0(2-l) and i3CO(2-l) emis- 
sion lines (see Table [TJ from the NGC 2264-C protocluster 
were taken with the IRAM-30 m telescope between October and 
November 2006 using the HEterodyne Receiver Array HERA 
( [Schuster et al.||2004| together with the VESPA autocorrela- 
tor backend. HERA is an assembly of two focal plane arrays 
(HERAl and HERA2) of 9 SIS receivers each. The two ar- 
rays have orthogonal polarization and their 9 pixels are arranged 
in the form of a center-filled square and are separated by 24". 
HERAl and HERA2 were used simultaneously, allowing us to 
map the same region of the sky at two different frequencies 
at the same time. VESPA was used with a channel spacing of 
320 kHz for 1^00(2-1), and 80kHz for C''^0(2-l) and '^CO(2- 
1). The VESPA autocorrelator configurations used during our 
observations can be found in Table [T] Observations were car- 
ried out in position switched On-The-Fly (OTF) mode, scanning 
NGC 2264-C in right ascension and declination. We adopted a 
drift speed of 2"/sec with a dump time of 2 s, and HERA was 
rotated by 9.5 ° to obtain a spacing of 4" between adjacent scan 
lines, while a derotator was used to keep the HERA pixel pat- 
tern stationary in the Nasmyth focal plane. All this resulted in a 
sampling of 4" over the entire map. The half power beamwidth 
(HPBW) of the 30-m telescope being of 11" at this frequency, 
this corresponds to an over-Nyquist sampling. The resulting map 
has a size of 3.3' x 3.3' (equivalent to ~ 0.8 pc^ at the distance 
of the protocluster). The emission-free reference position was 
located at offsets (-1500", -1500") from the protostellar core 
C-MM3, located at the center of our map. The typical system 
temperatures varied between ~300 K and ~600 K. The telescope 
pointing was checked every ~2 hours on J0528H-134, while fo- 
cus checks and adjustments were made every ~3 hours on the 



2.2. Data reduction 

All data were reduced with the CLASS 90 program, a new part 
of the GILDAS software package. Dedicated CLASS90 scripts 
were written to automatically reduce the data and produce a dat- 
acube for each of the three observed transitions (see Table[T]l. 
Linear (first-order) baselines were determined from velocities 
ranges where emission is no longer positively detected above the 
3cr noise limit anywhere in the map, and then subtracted from the 
spectra. We stress that our HERA observations resulted in very 
good baselines (see example in Fig. ffl. We used 5 km.s^'-wide 
velocity windows at either end of the CO(2-l) spectra to fit the 
substracted baselines, making us sensitive to any outflow emis- 
sion present up to LSR velocities of Vsyji 55 km.s"' . 
Some velocity channels showed spikes due to band splitting of 
the VESPA autocorrelator We had to flag them and replace their 
values by interpolating the closest two good channels. The worst 
case occurred around Vlsr ~ +23 km/s with four consecutive 
bad channels, leading to an interpolation over ~ 1.7 km/s, which 
we consider as being negligible for the present study of high- 
velocity emission. 

Finally, in order to improve the signal-to-noise ratio, we also 
smoothed the initial velocity resolution by a factor of two for 
the '^CO(2-l) map, leading to a final spectral resolution of 0.84 
km/s per channel. 



3. IVIapping results and analysis 

We detected a total of eleven sub-regions or "lobes" exhibiting 
high-velocity emission in the i2cO(2-l) map (Fig. ^. These 
eleven lobes are spatially distributed around the millimeter con- 
tinuum cores identified by Peretto et al. (2006 , 2007) , and four of 
these lobes can be directly associated with Class - like objects, 
as discussed in § 5.1 below. 
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Fig. 2. ^^CO(2-l) map of the NGC 2264-C protocluster. Blue contours show the levels of '^CO(2-l) intensity integrated between 
-27 km/s and 2 km/s in the blue-shifted part of the line, and go from 5 to 98 K.km.s"' . Red contours are levels of intensity integrated 
between 13 km/s and 34 km/s in the red-shifted part of the line, and go from 5 to 110 K.km.s"' . The eleven outflows discovered in 



the present survey are labeled by Fl to Fll . The background greyscale image is the 1 .2 mm dust continuum map obtained by Peretto 
et al. ( 2006|l with MAMBO on the IRAM-30m. White crosses and numbers refer to the millimeter continuum peaks identified by 
Peretto et al.](p006l[2007l) (see also Fig.fl). 



In §3.1 we explore the kinematical properties of these eleven 
outflow lobes, while §3.2 discusses the optical depth of the 
'^CO(2-l) transition in our map. 

3.1. Spatial extent of the outflow lobes and minimum 
associated velocity intervals. 

In order to determine the spatial extent of each outflow, we need 
to find out which pixels exhibit high-velocity emission due to 
outflows. We know from previous work that the NGC 2264- 
C protocluster exhibits a complex pattern of cloud velocities: 
the mean LSR velocity of the cloud was found to be 7.5+0.2 
km/s, but a pronounced velocity discontinuity was also found 
by |Peretto et al.| (|2006) near CMM3, and local systemic veloci- 
ties can depart from this value by up to 0.7 km/s, depending on 
position in the protocluster. Likewise, the intrinsic '-CO(2-l) 
linewidth varies with position as column density and turbulent 



motions vary through the protocluster. Therefore, to discrimi- 
nate emission due to high-velocity outflow motions from emis- 
sion due to an intrinsically wider '^CO(2-l) line from the cloud 
itself, and thus identify those pixels with outflow emission, we 
adopted the following conservative method. 

We first constructed '-CO(2-l) reference spectra for the 
local ambient medium surrounding the high-velocity lobes. In 
practice, we constructed each reference spectrum by averaging 
a large number of spectra (from 30 to 200 spectra depending on 
the subregion considered) located in the close neighborhood of 
each candidate outflow. In this way, we obtained one reference 
spectrum per high- velocity lobe for outflows F3 to Fll. We con- 
structed two reference spectra for each of Fl and F2, so as to 
ensure that small variations in the systemic velocity occurring 
between different sub-regions of these extended lobes be prop- 
erly accounted for In a second step, we compared each reference 
spectrum with the corresponding high-velocity lobe line profiles. 
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Fig. 3. Examples of high-velocity '^CO(2-l) spectra showing 
emission at blue-shifted (a) and red-shifted (b) velocities for 
the candidate outflows F8 and F9, respectively. In both cases, 
the dashed spectrum shows a '^CO(2-l) reference spectrum ob- 
tained by averaging ~100 spectra from the ambient cloud emis- 
sion in the immediate vicinity of the outflow. Vertical solid lines 
highlight the velocity range over which emission is no longer 
detected from the ambient cloud, whereas blue-shifted or red- 
shifted high-velocity emission is still present in the spectrum 
(above the 3cr noise level shown as an horizontal dot-dashed 
line) 



o 
o 
o 



9''28'40" 



1 1 
Outflow lobe 


1 

CMM1 1 

/ \ ■ 








1 \ 



15 8 



a (J2000) 



18:0 



4 














; Bullet ? 








\ : 







1 





-50 



Velocity (km/s) 



50 



Fig. 4. Spatial extent of the red-shifted outflow lobe F4. The 
upper panel shows the high-velocity '^CO(2-l) emission, in- 
tegrated from -1-14.3 km.s"' to -1-23 km.s"' (velocity interval 
marked by dashed lines in the lower panel). In this case, the 
reference spectrum associated to F4 does not show any signif- 
icant emission at LSR velocities larger than 14.3 km.s"^ The 
^cTspec level mentioned in the text corresponds to the thick black 
contour. All pixels showing integrated emission above this value 
were therefore assigned to outflow F4. A typical '^CO(2-l) 
spectrum observed near the peak of outflow emission is shown in 
the lower panel. Note the quality of the baseline and the detached 
high-velocity feature at Vlsr — 1-18 km.s"' . The area over which 
such a "bullet-like" emission feature (see § 5.1) is seen in F4 is 
marked by the white contour in the upper panel. 



Each pixel which showed emission stronger than three times the 
typical noise level, cTspec, in the '^CO(2-l) outflow spectrum in 
three consecutive channels (corresponding to a 2.5 km/s-wide 
window in velocity), at velocities where emission was no longer 
detected in the reference spectrum, was assigned to the candidate 
outflow under consideration. 

Such cases are illustrated in Fig. [3] using spectra from the 
two outflows F8 and F9, for which obvious high-velocity wings 
are observed, at blueshifted and redshifted velocities, respec- 
tively. The vertical solid lines mark the velocity range over 
which emission is only due to outflow material: in this velocity 
range no emission is detected in the reference spectrum, while 
emission is positively detected in the spectra from the candi- 
date outflow. Here, we ignore any additional outflow emission 
at lower velocities where emission is still seen in the reference 
spectrum. Distinguishing between low-velocity outflow emis- 
sion and turbulent broadening of the line spectrum from the am- 
bient cloud is more uncertain. We will try to account for the low- 
velocity outflow emission when we estimate outflow parameters 
in § 4 below. For the purpose of determining the spatial extents 
of the outflows, however, only those pixels satisfying the selec- 
tion criteria described above were included. 



Following this procedure for all pixels in the map, we were 
able to determine the spatial extent of each outflow lobe (see 
Col. 4 of Table |2]for the values, and Fig.|4]for an example). The 
velocity range over which high-velocity emission was detected 
from the outflow lobe but no emission was seen in the reference 
spectrum will be adopted as the main LSR velocities interval to 
be considered in further calculations. In the following, |Vi| will 
denote the extreme LSR velocity observed toward the outflow 
(beyond which To„,y7on, < lia-spec), and IV2I will denote the min- 
imum LSR velocity at which unambiguous outflow emission is 
detected (beyond which T,,,/ < lia-spec while T outflow > ^o-spec)- 
As can be seen in Table |2] the velocity ranges spanned by the 
observed high velocity emission vary from outflow to outflow, 
and '^CO(2-l) emission is detected at very high LSR velocities 
in some cases (up to -27 km/s in the case of the blue-shifted 
lobe F2). 

3.2. Optical depth of the ^^C0(2-1) emission 

In order to derive the masses and moments for the eleven can- 
didate outflows, we need to assess the optical depth of the CO 
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Table 2. Observed properties of the '^CO(2-l) outflows detected in the map of Fig. 2 



OUTFLOW 


Local systemic 


Outflow main 


Area 


Length 


LOBE 


velocity Vsy, 


LSR velocities interval 








(km.s-') 


(km.s 

Vi 


') 

Vz 


(arcsec^) 


(arcsec) 


Fl (Red) 


7.5 


28.5 


13 


2780 


127 


F2 (Blue) 


7.6 


-27 


1.1 


3060 


184 


F3 (Blue) 


7.4 


-26.6 


1.4 


1935 


70 


F4 (Red) 


7.7 


23 


14.3 


900 


53 


F5 (Blue) 


7.3 


-12 


0.3 


1080 


50 


F6 (Red) 


7.1 


28.1 


13.1 


1360 


56 


F7 (Red) 


7.2 


33.2 


14.7 


2230 


72 


F8 (Blue) 


7.2 


-23 


1.7 


590 


38 


F9 (Red) 


7.2 


30.2 


13.7 


800 


53 


FIO (Red) 


7.2 


32.2 


13.7 


610 


50 


Fl 1 (Red) 


7.7 


33 


15.2 


4180 


90 



Systemic LSR velocity of the ambient material around each outflow. 

Velocity interval [Vi ; V2] over which each outflow emission is detected but no significant emission is seen in the reference spectrum. 
Surface area of each outflow lobe. 
Major size of each outflow lobe. 



emission. We used our observations in three isotopes (see Table 
[T]) to estimate the optical depth of the CO(2-l) transition in our 
data. 

We first used our C''^0(2-l) data to check that the '^CO(2- 
1) emission observed away from line center was optically thin 
over the entire mapped area and estimated a maximum optical 
depth t('^CO) ~ 0.75 at the systemic velocity of the ambient 
cloud. 

We then derived the ^^CO(2-l) optical depth in the high- 
velocity gas from the observed '^CO/'^CO line intensity ratio 
(cf. Goldsmith et al.|1984| l. We measured this ratio at eleven po- 
sitions, chosen to coincide with the positions of the candidate 
outflows in the '^CO map. Assuming that the excitation temper- 
atures of the '^CO and '^CO lines are equal, and knowing that 
the '^'CO emission is optically thin, the '^CO optical depth '^t(v) 
at any given velocity v is related to the '^CO/'^CO intensity ratio 
by: 



1 - exp 



"t(v) 



1 - exp 



't(v) 



12t(v) 



(1) 



where X,-5o is the ['2CO]/['^CO] abundance ratio, here assumed 
to have the terrestrial value of 89 ( [Wilson & Matteucc"i||1992[ l. 
The signal-to-noise ratio in the '^CO(2-l) spectra was not ade- 
quate to derive '^Tv over the full velocity extent of the line wings 
seen in '^€0(2-1). We calculated the isotopic ratio up to ve- 
locities where the '^CO(2-l) emission became lower than the 
3cr noise limit, on either side of the line centroid for the eleven 
outflows. Figure ID shows the '2cO(2-l) to '^CO(2-l) ratio for 
three different outflows (outflow F4 as a solid line, outflow F7 as 
a dashed line, and outflow F3 as a dotted line). 
The '2CO(2-l)/'-'CO(2-l) intensity ratio is expected to increase, 
and '^Tv to decrease, at larger velocity offsets from the line cen- 
troid, but the '■'CO(2-l) emission becomes undetectable at these 
high velocities. At the highest velocities showing detectable 
'^CO(2-l) emission, a maximum '2cO(2-l)/'^CO(2-l) inten- 
sity ratio of ~ 45 is measured for both blueshifted and red- 
shifted gas (cf. Fig. [sjl, which corresponds to '^rn,,„js ~ 1.5. 
Therefore, '^t is much lower than 1.5 at LSR velocities larger 
than 12 km/s or smaller than 2 km/s. As the velocity intervals 
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Fig. 5. '2CO(2-l) to '^CO(2-l) intensity ratio calculated over 
the widest velocity range for which the '^CO(2-l) intensity is 
larger than the 3cr noise limit. Outflow F4 is represented as a 
solid line, outflow F7 as a dashed line, and outflow F3 as a dot- 
ted line. The vertical axis on the right shows the corresponding 
values of the '^00(2-1) optical depth, '^t(v). 



derived in § 3.1.2 for the candidate outflows lie outside the [2 : 
12] km/s range, we conclude that the '^00(2-1) outflow emis- 
sion is optically thin at both blue-shifted and red-shifted veloci- 
ties. 

Note, however, that we will apply an optical depth correc- 
tion factor when we derive the amount of mass entrained by the 
outflows at low velocities (i.e., inside the [2 : 12] km/s range) 
in § 4.1.2 below. For this purpose, we will use the values of 
'2t(v) estimated at LSR velocities between -1-2 and -1-12 km/s 
(see Figure[5]for examples). 



4. Momentum injection rate due to outflows in 
NGC 2264-C 

Our '^€0(2-1) map revealed the presence of eleven outflows in 
the NGC 2264-C protocluster. To quantify the effective feedback 
that these eleven outflows produce on the protocluster, we now 
estimate the masses of outflow entrained material, as well as the 
outflow dynamical parameters. 
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4.1. Outflow mass estimates 
4.1 .1 . Method of derivation 

Assuming local thermal equilibrium (LTE), the total gas mass of 
an outflow can be obtained from the observed '^CO(2-l) spectra 



as follows (see, e.g., |Scoville et al.| ( [T986j )): 



Te, + 0.93 



xd X 



Jv„ J A 



) dv.dA (Mo), 
(2) 

where T^x is the CO excitation temperature (in Kelvin), and d 
the distance to the source expressed in kpc. /(v) is related to 
the observed antenna temperature T^{\) (expressed in Kelvin) 

by /(V) = [j^] X T;(v), and j^^ /(v) d^.dA (with /(v) in 

K, dv in km.s"', and dA in arcsec ^) is the '^CO(2-l) intensity 
integrated over the surface area A, and over the velocity interval 
[Va : Vb] of the outflow. In deriving this formula, we adopted 
a standard CO to H2 abundance ratio, [^] = 10 (Frerking, 
Langer, & Wilson 1982). 

4.1 .2. Masses of outflow-entrained gas 

We first computed the minimum mass of gas associated with 
each outflow, assuming that the whole outflow emission is in- 
cluded in the main LSR velocity range [Vi : V2] derived in 
§ 3.1.2 (Col. 3 of Table [2|. In this first estimate, we thus did 
not take into account any additional outflow emission seen at 
lower absolute velocity offsets from line center, where emission 
is positively detected in the reference spectrum (see Fig.|6]l. We 
showed in § 3.2 that the '^CO(2-l) emission is optically thin ev- 
erywhere outside the velocity range [2 : 12] km/s in our map. As 
the minimum velocity intervals considered in this first estimate 
all lie outside the [2 : 12] km/s range, the minimum masses of 
outflowing gas can thus be computed as: 



1.6 X 10"" X 



/7 

JVi J A 



Tliw) dw.dA (Mo) (3) 



In the above equation, r*(v) is expressed in Kelvin, d\ is a 
velocity bin (in km.s"'), while dA corresponds to the surface 
element (in arcsec^). 

The main source of uncertainty on the derived outflow-entrained 
masses is the uncertain value of the CO excitation temperature. 
For optically thin CO emission, the excitation temperature 
is expected to be close to the gas kinetic temperature due to 
the low dipole moment of the CO molecule. Nevertheless, 
large temperatures (up to 80 K), due to gas heating by outflow 
shocks (Umemot o et"aL][T992l |Bachiller et aL][T99^ |2001i ) are 
sometimes found in CO outflows. As can be seen from Eq. (2), 
the derived outflow mass is almost proportional to the assumed 
excitation temperature. We assumed an excitation temperature 
of 20 K in our calculations, which is representative of typical 
CO ouflows. The uncertainty in the CO excitation temperature 
introduces a factor of ~ 2 uncertainty on our estimates of the 
masses, momenta, and momentum fluxes. 
We calculated the value of " ^I*-^^ dv.dA directly from our 
data by integrating the intensity of the '^CO(2-l) emission, 
both spatially over the outflow extent, and kinematic ally over 
the main LSR velocities range [Vi ; V2] of each outflow. The 
outflow extents and main LSR velocities ranges adopted are 
given in Table |2] 




-10 , , , 

Velocity (km/s) 

Fig. 6. Illustration of the velocity ranges used for the calculation 
of minimum and maximum masses of the outflows. The dashed 
spectrum shows the '^CO(2-l) reference spectrum from the am- 
bient cloud. The solid spectrum corresponds to the candidate 
outflow F3. Vertical solid lines show the velocity ranges used in 
the computation of the outflow masses. [Vi ; V2] represents the 
interval used for the calculation of the minimum outflow mass 
(see §4.1.2 for details). [V2 ; V3] is the additional velocity inter- 
val used to compute Madd by integrating the excess '^CO(2-l) 
emission (area highlighted in grey in this case). 



The minimum mass of gas associated with each outflow is given 
in the fifth column (first value) of Table |4] The minimum to- 
tal mass of gas entrained by the network of eleven outflows in 
NGC 2264-C is estimated to be 0.3 + 0.15 Mo. 

In order to estimate the maximum mass of gas entrained by 
the eleven outflows observed, two calculations were made. 
First, for each outflow lobe, we assumed that all of the blue- 



outflow 



>Tref + 30-. 



spec 



shifted or red-shifted emission such as T, 
was due to the corresponding outflow (see Fig.[6]l. This hypoth- 
esis may overestimate outflow masses by including low-velocity 
motions due to local turbulent broadening of the ambient cloud 
emission. The corresponding values therefore provide upper lim- 
its to the masses of outflow-entrained gas. We also recall that, 
at low velocities, the '^CO(2-l) emission is no longer optically 
thin (see Fig. [5]). Therefore, we used the values of optical depth 
derived previously to correct for optical depth effects on our 
mass estimates. The additional outflow emission at low veloc- 
ities (see Col. 3 of the online Table [3] for numerical values) can 
be expressed as: 



Madd = i.exio^v 



r r— 

Jv, J A 1 - e---' 



[n(v)|flo„-Tl(v)|,ef]dv.dA (Mo), 



(4) 



3cr, 



spec 



where IV2I is the velocity at which r^(v)|ref < 3crref 
and IVgl the velocity at which r;(v)|flow = r;(v)|rf - 
(if T'KV.yJIflow > T;(V,y,)|,ef then IV3I = V,,,). 
As a second step, we further took into account the material en- 
trained by outflows down to velocities comparable to the sys- 
temic velocity of the surrounding medium. We first extracted a 
characteristic mean spectrum per outflow, and fitted a power law 
to this mean spectrum, Toutflow - /3 x (v - VsysY', in the veloc- 
ity range Vi < v < V3. We then extrapolated the fitted profile 
for velocities V3 < v < Vjys, in order to estimate the amount of 
mass entrained at very low velocities and which remains hidden 
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in the main body of the '^CO(2-l) line: 



§3.1). This results in an inclination correction factor 



Mio„ = 1.6xlO"V 



Jv3 J A 1 - 



(5) 

where Vsys is the systemic velocity associated with the outflow 
(see Col. 2 of Table[2|. 

The power-law indices of the fit can be found in the Col. 4 of 
the online Table[3] The additional mass of outflow-entrained gas 
at velocities between V3 and the systemic velocity of the cloud 
Vsys is given in Col. 5 of the online Table [3] We obtained satis- 
factory fits for all eleven outflow lobes, except outflow F4 which 
exhibits a detached high-velocity emission feature (see Fig. |4]i. 
The best fits lead to a fraction of invisible, low-velocity outflow- 
entrained mass Miow/A^max ~ 10% - 45%. Since this additional 
mass is moving at very low velocities (1 to 3 km.s"' away from 
the systemic velocity), its contribution to the outflow dynamical 
parameters will nevertheless be quite small (see § 4.2 below). 

Finally, the maximum masses of outflow-entrained gas M^ax 
were computed by adding the masses M-^^d + -^low calculated in 
the velocity interval [V2 ; Vsys] to the minimum masses M,„in 
derived earlier. These maximum masses are reported in Col. 5 of 
Table |4] (second value). We stress that these masses represent the 
maximum possible values of the outflow masses, as the hypothe- 
sis made in their derivation led us to include all of the '^CO emis- 
sion potentially due to outflows. The overall uncertainty on each 
of the outflow mass components (Mj^jn, M.^^^, Mio^v) is about a 
factor of two. 



4.2. Outflow dynamical parameters 

Using the mass estimates given in Col. 5 of Table [4j we first 
computed the momentum associated with each outflow: 



flow 



= Mflow X V 



char 



(6) 



where Vchai is the characteristic velocity of the outflow lobe, 
computed as the intensity-weighted mean absolute velocity of 
the outflow over the relevant velocity range : 



,/obs 



Vchar - 



1 



cos(() A',„,-Xcos(r) 



pixels 



2 T*(v)|v-V,y,| 

Av 

E T*(v) ■ 



The inclination angle / between the flow axis and the line of 
sight (1.0. s.) has to be taken into account when computing the 
true characteristic velocity Vchar, since it is related to the ob- 
served characteristic velocity (projected on the l.o.s.) V*^^. (see 
Col. 2. of Tableg, by Vchar = V°^^^ / cos(i). 
We derived the outflow momentum in each of the three velocity 
regimes [Vi ; V2], [V2 ; V3], [V3 ; Vsys] considered for the cal- 
culations of the outflow -entrained mass (cf. § 4.1.2 and Fig. |6]l. 
The total outflow momentum was calculated as the sum of the 
three momenta estimated in the three velocity regimes. 

The outflow momentum flux was obtained by dividing the 
momentum by the estimated dynamical timescale t^yn of the out- 
flow: 



-^flow ~ -f*flow/^dy 



(7) 



The dynamical timescale, fdyn = ^flow/Vmax, also depends on 
the inclination angle /. Indeed, / affects both the estimate of the 
outflow length (/flow = / / sin(/) where I is the observed pro- 
jected length on the sky), and the velocity estimates since the 



highest intrinsic flow velocity is Vmax^Vj^^^ / cos(/), where 
Vjjgx is the highest observed velocity (corresponding to Vi in 



cos(0 
sin(0 



for 



the dynamical timescale. Given the inclination correction on the 
momentum, the global inclination cor rection factor for the m o- 

mentum flux is thus f(i) = '^'"ff, (cf. 

■' ^ ' cos^(0 ^ 



Bontemps et al. 



19961. 



The results of these calculations are given in Table |4] The 
length of each outflow is reported in Col. 3, where the first 
quoted value is the observed projected length of the outflow 
and the second value is corrected for the mean inclination angle 
i= 57.3°. Column 4 gives the dynamical timescale of each out- 
flow. The lower value is the apparent dynamical timescale (with 
no correction for inclination), while the upper value is corrected 
for the mean inclination angle i= 57.3°. The outflow momentum 
and momentum flux values are given in Col. 6 and Col. 7. The 
minimum quoted values were derived using the minimum mass 
estimates and ignoring inclination effects, i.e. assuming that the 
outflows had their axes along the l.o.s. The maximum values 
were derived using the maximum mass estimates and assuming 
a mean inclination angle i~ 57.3° to account for inclination ef- 
fects (leading to a mean correction factor {/(/)) = 2.9 for the 
momentum flux - cf. |Bontemps et al.|1996| l. 

We stress that the amount of outflow mass hidden at low ve- 
locities in the main body of the '^CO(2-l) line, Miovi,, contributes 
only ~ 1% to 5% of the maximum momentum flux computed 
here, due to the low velocity of this material. Therefore, thanks 
to the procedure adopted here to take into account the different 
velocity regimes of the outflows, the relatively large uncertainty 
in the outflow-entrained mass at low velocities results only in a 
small uncertainty factor for the momentum flux. 

The main source of uncertainty in our momentum and mo- 
mentum flux estimates comes from the hypothesis made on the 
inclination angle with respect to the l.o.s. . We considered that 
our sample of eleven outflow lobes was large enough to justify 
a statistical treatment and applied a mean inclination angle of 
i~57.3° assuming random outflow orientations. We estimate that 
this hypothesis on the inclination angle leads to a factor of ~ 2 
uncertainty on the individual momentum and momentum flux 
values. Altogether the global uncertainty on the derived momen- 
tum fluxes is roughly a factor of ~ 8, which is consistent with 
the results of [Cabrit & Bertout ( 1990p based on an analysis of 
synthetic outflow data. 

Table |5] gives the total dynamical parameters of the entire 
population of outflows observed in NGC 2264-C obtained by 
summing the values computed for the eleven outflow lobes (from 
Table |4]l. The total kinetic energy contained in the eleven out- 
flows was estimated from the total momentum as follows: Etot= 



X V, 



char 



3.4x10^4 - 1x10*5 ergs. We also es- 



2: 1/2 X Pflo 

flows 

timated the rate of kinetic energy deposited by the outflows in 
the surrounding medium as: Ltot= 1/2 x Ftot x Vesc, where Vesc 
is the e scape velocity of the clou d Vesc= V2GMcioud/^cioud ~ 5 
km.s ' (Stojimirovic et al. 2006 1. In this way, we found a total 



mechanical luminosity Ltot ~ 0.21 - 1.2 L©. 
The range of values we derive for the total mass of outflow- 
entrained gas is in agreement with the mass values found by 
Margulis et al. (1988) for the whole NGC 2264-C region. The 
range of values we find for the overall momentum and mechan- 
ical luminosity are in agreement only with the lower values of 
|Margulis et al. ( 1988[ l. This difference is mainly due to their 
method, which used the maximum projected flow velocity to 
compute upper limits to the momentum and mechanical lumi- 
nosity. Furthermore, Margulis et al. (1988i did not resolve the 
network of outflows in NGC 2264-C and thus used a single char- 
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OUTFLOW 


Characteristic velocity 


Length 


Dynamical time 


Mass 


Momentum 


Momentum flux 


LOBE 


"wobs 
char 


/ = 90° - i = 57.3° 


= 90° - i = 57.3° 










(km s"h 


VP*-/' 






CM^ km s^h 


dO^^ km s"' vr"'") 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Fl (Red) 


+9.5 


0.49 - 0.58 


24.5 - 15.9 


3.7-35.7 


0.58- 1.23 


2.4- 10.2 


F2 (Blue) 


-13 


0.71-0.84 


20.0- 13.0 


2.9 - 36.8 


0.37 - 0.97 


1.8-6.8 


F3 (Blue) 


-11.5 


0.27 - 0.32 


7.7 - 4.9 


4.2 - 34.2 


0.49 - 1.07 


6.3 - 28.3 


F4 (Red) 


+7.6 


0.2 - 0.24 


12.8 - 8.2 


0.8-3.1 


0.06-0.14 


0.5 - 3.2 


F5 (Blue) 


-12 


0.18 - 0.21 


9.1 - 5.8 


0.9 - 37.9 


0.12 - 0.74 


1.3 - 11.5 


FO (^KeQ^ 


+0.0 


U.Z - U.Z+ 


y.j - J.y 


i .0 - J J.J 


U. i - U. / 1 


11 1 n ^ 


F7 (Red) 


+10.3 


0.28-0.33 


10.6 - 6.8 


6.0- 111 


0.62 - 2.56 


5.8-35.3 


F8 (Blue) 


+10.9 


0.14-0.16 


4.6 - 2.9 


0.5-7.8 


0.05-0.16 


1.2-4.4 


F9 (Red) 


+11 


0.2 - 0.24 


8.8-5.6 


0.6-6.8 


0.06-0.18 


0.7 - 3.4 


FIO (Red) 


+11 


0.18-0.21 


7.9 - 5.0 


0.5 - 7.4 


0.05-0.18 


0.7 - 3.4 


Fll (Red) 


+ 10.8 


0.35-0.41 


13.5 - 8.64 


11.2-93.2 


1.21 - 3.12 


8.9-55.3 


L1551(" 




0.65 - 0.77 


48-31 


530 - 670 


18-23 


47 - 136 


LI 157 (i~81°)<" 




0.4 


2.5 


62 


31 


1100 



; Vsys] velocity range 



Observed characteristic offset velocity of the outflow lobe (see details in §4.2), computed over the [V 
Length of each outflow lobe. 
Dynamical timescale of each outflow. 

Minimum and maximum masses of outflow-entrained gas, computed as explained in § 4.1. 
Minimum and maximum momenta, computed as explained in § 4.2. 
Range of momentum flux derived for each outflow (see § 4.2. 

The last two rows give, for the sake of comparison, the values computed for the two bipolar outflows LI 55 1 and LI 157, used by 
[&Lilp007l l. 



Nakamura 



Table 5. Global properties of the population of eleven '^CO(2-l) outflows observed in NGC 2264-C. 



M,„, 


Plot 


Fiot 


Etot 


Ltot 


(10-^ Mo) 


(Mo.km.s-') 


(10^^ Mo.km.s-'.yr-') 


(10'*'' ergs) 


(Lo) 


(1) 


(2) 


(3) 


(4) 


(5) 


3-189 


3.8-11 


30-170 


3.4-10 


0.2-1.2 



' ' Minimum and maximum total mass of outflow-entrained gas, computed from the individual outflow masses given in Tabled 

Minimum and maximum total outflow momentum, computed from the individual outflow momenta given in Tablepl(Pto, = 2 Pfiow)- 

M flows 

Minimum and maximum total outflow momentum flux (F,o, = 2 Ffiow)- 

flows 

Minimum and maximum total kinetic energy of the outflows. 

Minimum and maximum total mechanical luminosity of the outflow population. 



acteristic velocity for the eleven outflows, while our study re- 
veals large velocity variations from outflow to outflow. 



5. Discussion 

5.1. Identification of driving sources 

The lobes named Fl and F2 in Fig. |2] are the most elon- 
gated and best coUimated outflow lobes identified in our study. 
Fl extends over 130" and F2 extends over 180", corresponding 
to physical lengths of 0.49 pc and 0.7 1 pc, respectively, assuming 
the outflows are perpendicular to the line of sight (i= 90°). The 
similar morphologies and relative positions of the two lobes in 
the map of Fig. |2] strongly suggest that Fl and F2 are physically 
associated and correspond to the two lobes of a single bipolar 



outflow. Based on Fig. [2] and Fig. [7] the best candidate driving 
sources for this bipolar flow are C-MM3 and C-MM13. 

To determine which source is driving the (Fl, F2) bipolar 
outflow, we constructed several position-velocity (PV) diagrams 
along the axis of the flow, assuming various barycenters and ori- 
entations for the outflow. Figure|8]shows the best (most symmet- 
ric) position-velocity diagram we found. This PV diagram was 
constructed along an outflow axis passing through C-MM13 and 
defined by a position angle of 12° (east of north) on the plane of 
the sky. Note that the PV diagrams constructed for position an- 
gles between 12° and 22° are nearly indistinguishable. Figure |7] 
indicates the range of possible outflow axes based on our PV- 
diagram study. It is clear from Fig. |8] that the highest- velocity 
emission (+20 km/s with respect to the systemic velocity) is 
detected in the immediate vicinity of the protostellar core C- 
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Fig. 7. Blow up of the '^00(2-1) map of Fig.|2]toward the (Fl, 
F2) bipolar outflow in the central part of the NGC 2264-C pro- 
tocluster. The full contours show the integrated intensity in the 
blue-shifted part of CO(2-l) line (from 1 to -27 km/s). Contours 
range from 5 to 98 K.km.s The dashed contours are levels 
of CO(2-l) integrated intensity in the red-shifted part of the 
line (from 13 to 28.5 km/s), and contours range from 5 to 110 
K.km.s"'. Markers refer to the positions of millimeter contin- 
uum peaks C-MM3 and C-MM13. The solid line arrow marks 
the direction of the bipolar outflow made up by Fl and F2, at 
a position angle of 12° (east of north). The dashed line arrow 
marks the direction of the (Fl, F2) bipolar outflow, at a posi- 
tion angle of 22°. This corresponds to the range of outflow axes 
for which the most symmetric PV diagrams are obtained (see 
Fig.lHl. 



Flow with P.A. of 12° : 



20 




-20 



Levels (K) : 
0.2 0.4 0.8 1.4 2 



-100 100 

Offset from C-MM13 along outflow axis (") 

Fig. 8. Position- velocity diagram of the (Fl, F2) bipolar flow, 
constructed along an axis passing through C-MM13 (offset 0") 
and defined by a position angle of 12° in the plane of the sky, and 
oriented from south to north (see Fig. |7]i. The dotted horizontal 
line represents the value of the systemic velocity observed at the 
center of the protocluster. The dashed line in red represents the 
best symmetry axis for the highest velocity emission. Note that 
most of the red-shifted velocity occurring at negative offsets (i.e. 
southern) is due to outflow Fll and is unrelated to (F1,F2). If 
we ignore this Fll emission feature, the red-shifted and blue- 
shifted high-velocity components seen in the diagram exhibit a 
symmetric pattern with respect to offset 0" (intersection of the 
dotted and dashed lines), which corresponds to the position of 
the protostellar core C-MM13. 



MM 13, and that this high-velocity emission presents a very high 
level of symmetry about C-MM13. On this basis, we conclude 
that C-MM13 is the most likely driving source of the bipolar 
outflow made up by lobes Fl and F2. 

One possible explanation for the slight misalignment be- 
tween the outflow features located close to C-MM13 and their 
more distant counterparts (see Fig. |7]i is that the protostellar 
jet driving the (Fl, F2) outflow precesses with time (Gueth et 



al. 19981). High-velocity features located farther away from the 



driving source, which presumably correspond to older ejection 
events, would then indicate a different outflow direction than 
that indicated by the high-velocity features located closer to C- 
MM13. 

Based on spatial coincidence and absence of confusion (see 
Fig. |2|, we also believe that the blue-shifted lobe F3 is driven 
by the millimeter source C-MM2 and that the red-shifted lobe 
F4 is related to C-MMl. The seven other lobes are either too 
poorly coUimated (eg. red-shifted lobe Fll), too far away from 
the protostellar clump (eg. red-shifted lobe F6), or too blended 
with other lobes to be directly associated with any particular 
Class 0-like object. 

The outflow lobe F4 exhibits a compact (~ 0.03 pc) '^CO 
emission peak at Vlsr — 1-18 km.s"', corresponding to a pro- 
jected velocity offset — nil km.s"' from the systemic velocity 



(see Fig. |4|l. No emission is seen at Wlsr - -1-18 km.s"' in our 
'^CO data, ruling out a background cloud origin for the '^CO 
peak. Such a discrete '^CO emission feature, well defined in both 
space and velocity, is reminiscent of the high-velocity molecu- 
lar "bullets" observed toward some young outflows driven by 
Class objects such as the L1448-C bipolar flow (Bachiller et 
|aL] [T990l ). 

We also stress that there is evidence for at least two well- 
coflimated outflows in NGC 2264-C: the first one is the (Fl, F2) 
pair which is very likely driven by CMM-13, while the second 
one is the outflow lobe FIO which is also well collimated but of 
unidentified origin. Given that the best collimated CO outflows 
are generally detected toward the youngest protostellar objects 
( [Andre et al.|1990||Bachiller|1996 ), these highly-collimated out- 
flows further testify to the presence of extremely young (Class 0) 
protostellar objects in the NGC 2264-C protocluster. 



5.2. Properties of the driving sources and comparison witii 
well-studied outflows 

We now discuss the properties of the three protostellar driving 
sources identified in the previous section : C-MMl, C-MM2, and 
C-MM13. 
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5.2.1. Source properties 



We used 5/?/fzer/MIPS data at 70 //m, APEX/P-ArTeMiS data 
at 450 fi m (see Fig. \\\)), an d IRAM 30m MAMBO data at 
1.2 mm ( jPeretto et al.| pOO F) to estimate the bolometric lumi- 
nosities (Lboi) of C-MMl, C-MM2 and C-MM13. We could di- 
rectly measure 70 fim and 450 fim flux densities for the rela- 
tively isolated sources C-MMl and C-MM2 and we then esti- 
mated their bolometric luminosities using the online tool devised 
by Robitaille et al. (2007 ) to fit the spectral energy distributions 
of young stellar objects (YSOs). Unfortunately, due to the pres- 
ence of the bright young star IRSl (see Fig.[T]i in the field, it was 
impossible to derive a flux from the Spitzer/MlPS data for C- 
MM13. We could nevertheless obtain a MIPS 70 i-im flux density 
and a bolometric luminosity for the neighboring source C-MM3 
(~13.5" away from C-MM13). We then derived a luminosity es- 
timate for C-MMl 3 assuming that the bolometric luminosity ra- 
tio L^^",*^'^'^/L^^",*^*^^ was equal to the 450 /im flux density ratio 
^450^m^^/^5i)^m^ ■ strcss that the bolometric luminosity de- 
rived in this way for C-MMl 3 is more uncertain than the Lboi 
estimates obtained for C-MMl and C-MM2. 

Our bolometric luminosity estimates for the three identified 
driving sources are listed in Table [6] We also give the enve- 
lope masses (Mgnv) of these protostellar objects as measured by 
[Peretto etaT] ( |2()07 1) with the IRAM Plateau de Bure interferom- 
eter on -3000 AU scales. 



Table 6. Properties of the sources : C-MMl, C-MM2, C-MM13, 
and C-MM3. 



The protostellar sources C-MMl, C-MM2, C-MM3 and C- 
MM13 were then placed in an Menv-i'boi diagram (see Fig. |9]l 
similar to the ones presented by Andre et al. (2000 2008| l. Model 
evolutionary tracks are shown, computed by assuming each pro- 
tostar forms from a bounded condensation of finite initial mass 
Menv(O) and has Lboi = GM^(t)M^^^(t)/R^(t) + L*(f), where 
/?* is the protostellar radius and the interior stellar lumi- 
nosity ( Stabler 1988 Hosowaka & Omukai 2008 1. The mass 
accretion rate and the envelope mass were assumed to be re- 
lated by Macc(0 - eMaw(t)/T, where e = 50% is the typ- 
ical star formation efficiency for individual cores (Matzner & 
McKee| [2000 ), and r - IQ^ yr is the characteristic timescale of 
protostellar evolution, leading to Macc(t) and Menv(0 functions 



declining exponentially with time (see [Bontemps et aT] 1996 1. 
The positions of C-MMl, C-MM2, C-MM3 and C-MMl 3 in 
this diagram demonstrate that they are bona fide Class ob- 
jects. Moreover, comparison with the model evolutionary tracks 
suggests that these objects will form stars with main-sequence 




CLASS I 
(Me„v < M.) 
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Fig. 9. Envelope mass versus bolometric luminosity dia- 
gram comparing the locations of the protostellar sources of 
NGC 2264-C (open triangles with error bars) with the positions 
of low-mass Class I (filled circles) and Class objects (open 
circles) (from Andre et al. 2000). Model evolutionary tracks, 
computed for various final stellar masses m{'""', are shown. The 
dashed and dotted lines are two Mf,„,, - Lboi relations marking the 
conceptual border zone between the Class (Menv > M*/e) and 
the Class I (Mgnv < M,^/e) stages, where e ~ 50% is the local 
star formation efficiency (see Andre et al.|2008 for details). 
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5.2.2. Comparison witli oilier outflows from the literature 

Here, we compare the properties of the three outflows (Fl, F2), 
F3, and F4 with those of several protostellar outflows from the 
literature. 

[Bontemps et al. (T996| carried out a study of outflow ener- 
getics for a sample of 45 low-luminosity YSOs including both 
Class I and Class objects. They found that outflows from 
younger sources were more powerful. The outflow momentum 
flux values we find here, namely a few 10"^ Mo.km.s^'.yr"', 
are very similar to the momentum fluxes of the Class outflows 
in the Bontemps et al. sample (see Fig. 10 1. The left panel of 



Fig. [To] plots outflow momentum flux versus driving source lu 
minosity, while the right panel shows Fflo„c/L versus Mgnv/Li 



bol 



0.6 



(see Bontem ps et al.| ( |19 96 ) for details). Figure 10 shows that C- 
MMl, C-MM2, and C-MMl 3 drive outflows which are typical 
of low- to intermediate-mass Class objects and closely follow 



the trend established by Bontemps et al. ( 1996 1. 



We can also compare the properties of the outflows mapped 
in NGC 2264-C with the two molecular outflows used as refer- 



ences by Nakamura & Li ( 2007 1 to estimate the typical outflow 



momentum to be injected in their numerical simulations. 
The first one is the well-documented outflow driven by the 
border-line Class I/Class object L1551 IRS 5. jStojimirovic 
|et al.| ( |2006| showed that the LI 551 IRS 5 bipolar outflow ex- 
tends over 1.3 pc and drives a total momentum of ~ 17.8 - 23 
Mo.km.s The dynamical parameters of the L1551 outflow, as 



masses M: 



final 



2-8 Mp 



re-estimated by us from the results of Stojimirovic et al. (2006) 
for consistency with the present study, are reported in Table [4] 
(see also Fig. 10 1. We note that the momentum flux injected by 
the L1551 outflow in the surrounding medium is of the same or- 
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Fig. 10. CO momentum flux Ffiow versus bolometric luminosity (left), and FgowC/L versus Menv/Lboi°'^ (right). The three candidate 
Class objects identified as the driving sources of the outflows (Fl, F2), F3 and F4 in NGC 2264-C are shown as blue stars. For 
comparison, the Class and Class I objects studied by Bontemps et al. (|1996 1 are shown as open and filled circles, respectively. The 
driving sources of the two molecular outflows taken as references by Nakamura & Li ( 2007| l are shown by black error bars. 



der as the momentum flux of the strongest outflow observed in 

the NGC 2264-C protocluster (F 11 ) . 

The second reference outflow used by [Nakamura & Li ( 2007| l is 
the one driven by the low-luminosity Class protostar L1157- 
MM (e.g. [Umemoto et al.| 



1992, BachiUer et al. ,2001). The 



dynamical parameters of the LI 157 outflow are also reported 
in Table |4] and used in Fig. [TO] It can be seen in Fig. [TO] that 
the LI 157 outflow is characterized by a very high momentum 
flux. None of the outflows mapped in the NGC 2264-C pro- 
tocluster drives that much momentum flux in the surrounding 
medium, and even the combined eff'ects of the eleven outflows 
in NGC 2264-C inject at most Ftot ~ 1.7x10"^ Mo.km.s^'.yr"' 
(see Table [5]), i.e one tenth of what is found for the LI 157 out- 
flow alone. Therefore, we conclude that LI 157-MM represents 
an extreme case of a Class outflow. 

We also note that the total momentum flux derived here 
(Table |5]) for the network of eleven outflows in NGC 2264-C 
is similar to the total outflow momentum flux found by Knee & 
Sandell] ( [20001 ) in the NGC 1333 protocluster. 



5.3. Outflows as a source of turbulence and support 

5.3.1 . Observed turbulence vs outflow contribution 

We now discuss whether or not the eleven outflows mapped in 
NGC 2264-C are likely to contribute a significant fraction of the 
turbulence observed in the protocluster 

From our '^CO mapping, we determined a one-dimensional 
velocity dispersion o-y ~ 1.7 km.s"', measured on the mean 
'^CO(2-l) spectrum over the NGC 2264-C clump (radius R ~ 
0.7 pc). If we place NGC 2264-C in the linewidth-size dia- 

[T98T) , 



gram representing Larson's first scaling law (Larson 



we find that NGC 2264-C exhibits a velocity dispersion about 
twice as large as what is expected from the classical linewidth- 
size relationship observed in giant molecular clouds (GMCs) 
( [Solomon etlT] [1987] ). This suggests that NGC 2264-C is 



more turbulent than the CO clouds considered by [Solomon et[ 
al. (1987). However, 'Heyer et al. (2008| recently re-examined 
Solomon's results with new ^^CO data and found the relation- 
ship cTy = y/izCE/S X R^^^, which shows that the normaliza- 
tion of Larson's first scaling law depends on the mean sur- 
face density S of the clouds. From the 1.2 mm dust continuum 
observations of [Peretto et al.[ P006[ l we estimate a mean sur- 
face density of ~ 3000 Mo.pc^~for the NGC 2264-C cluster- 
forming clump, which is much higher than the mean surface 
density ~ 170 Mq pc"^ of the CO clouds following the classi- 
cal linewidth-size relation. The location of NGC 2264-C in the 
cTy I R"^ - f(E) diagram of Heyer et al. (2008]) is nevertheless 
similar to that of the densest clou ds discussed by the se authors. 
We also note that both the N2H+ (Peretto et al.| 12006 ) and '^CO 



(this work) linewidths of NGC 2264-C are consistent with the 



log(Av) = 0.23 -I- 0.21xlog(R) correlation found by Caselh & 



Myers 



1995 ) for massive cloud cores in Orion. Therefore, al- 



though NGC 2264-C exhibits broader linewidths than those ex- 
pected from Larson's classical linewidth-size relation, it exhibits 
the same level of turbulence as that generally found in massive 
clouds and cores. 

Next, we can estimate the rate of turbulent energy dissi- 
pation in the NGC 2264-C clump using equation (7) of Macl 
"Low ( 19 99). From the 1.2 mm dust continuum map obtained by 
iPeretto et al. ( 2006| ), a total gas mass of 2300 M© can be derived 
within a radius of 0.7 pc. Using the velocity dispersion derived 
from our '^CO observations, we find : 



l/2xMcr2 



1.2 L„. 



(8) 



For comparison, the total mechanical power due to the eleven 
outflows in NGC 2264-C is Ltot = 0.7 + 0.5 Lq (see § 4.2). We 
conclude that the network of protostellar outflows observed in 
NGC 2264-C can contribute a very significant, if not dominant, 
fraction of the turbulence observed in the protocluster 
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5.3.2. Gravitational collapse vs outflow support 

In order to determine whether the force exerted by the outflows 
on NGC 2264-C can slow down or even stop the global collapse 
of the protocluster, we first estimate the force needed to keep the 
whole protocluster in hydrostatic equilibrium. 
The pressure gradient needed to balance the gravitational accel- 
eration agrav('") = GM(r)/r^ at radius r and keep a spherical 
clump with a mass distribution M(r) in hydrostatic equilibrium 
is: 



dP„ 



dr 



-G 



p(r).M(r) 



\2 



(9) 



this 



Assuming a density distribution of the form p(r) 

corresponds to a supporting pressure PgnviR) = GxM{R)^ /^nR'^ 
for the spherical shell of radius R enclosing a mass 
M(R) = Icr^RIG. 

Therefore, the total force needed to balance gravity at radius R 
in file NGC 2264-C protocluster is: 



Fgrav(«) = Pgrav(«) X AuR^ = G X M{Ry jlR^ 

The corresponding gravitational binding energy is 
E„.AR) = -G X M{RflR 



(10) 



(11) 



If we adopt Ry^QA pc and M{Ri) = 1600 Mq, 
"■^ '^".Q.km.s^'.yr"', 
50km.s-i ^"'"-1 - "^f^"- ^"'"-2 



we find : 



Fgrav =40X10 

and a„r 



.Myr"' ~ SOpc.Myr" 



If we adopt /?2= 0.7 pc and M{R2) = 2300 Mq, we find : 



Fgrav = 23x10"^ Mo.km.s-i.yr-1, 
and flgrav = 20 km.s^'.Myr"' ~ 20 pc.Myr"^. 



Comparing Fg^av to our upper estimate of the force ex- 
erted by outflows on the NGC 2264-C protocluster, Ftot - 
1.7xlO"^M0.km.s"'.yr"' (see §4.2), we find that the eleven ob- 
served outflows fall short by a factor ~ 25 to provide significant 
support against collapse in NGC 2264-C. 

We also note that the turbulent acceleration (cf. |Matzner| 
[2007] l : 



(12) 



is only aturb ~ 4.7 km.s^'.Myr"' at radius R = 0.7 pc, which is 
~ 4-5 times lower than the gravitational acceleration derived at 
the same radius. Therefore, the enhanced level of turbulence ob- 
served in NGC 2264-C (cf. § 5.3.1) is also insufficient to support 
the clump against global collapse. 

From an energetic point of view, the total kinetic energy (see 
Etot in Table [5]l derived in §4.2 for the eleven outflows is nearly 
three orders of magnitude lower than the gravitational bind- 
ing energy of the NGC 2264-C cluster-forming clump (Egrav ~ 
5.5x10'*'' ergs). Stojimirovic et al. (2006) showed that the ki- 



netic energy injected by the L1551-IRS5 outflow in the sur- 
rounding medium was comparable to the binding energy of the 
LI 551 cloud itself. The kinetic energy they derived from the 
L1551 outflow (1.5x10'*^ ergs) is comparable to our upper es- 
timate of the total kinetic energy of the outflow population in 
NGC 2264-C (see Table |5|, but the L1551 cloud is much less 
massive (110 Mq) than NGC 2264-C for an equivalent vol- 
ume (radius of 0.5 pc). Likewise, we find that the total mass of 
outflow-entrained material represents only 0.25% of the mass of 
the NGC 2264-C protocluster, whereas the IRS5 outflow repre- 
sents ~5% of the total mass of the L1551 cloud. 



For comparison, the rate of gravitational energy release due 
to global collapse is : 



dE„ 



dt 



= -(GM(R)2/R')r=|Eg,,, 



Vi 



infall 



R 



15 Lo, 
(13) 

in the NGC 2264-C, where we have used the large-scale infall 
velocity Vinfaii = 1.3 km/s measured by [Peretto et al. ( 2006| l. 
This is a factor of ~ 15-20 larger than both the total mechanical 
power of the observed outflows (cf. § 4.2) and the rate of 
turbulent energy dissipation in the clump (cf. § 5.3.1). 

We conclude that neither the direct force exerted by the net- 
work of observed protostellar outflows, nor the turbulence they 
produce, can support the cluster-forming clump against global 
collapse at the present time. Therefore, the additional source of 
support against gravity, not included in the purely hydrodynamic 
simulations of [Peretto et al.| ( |2007| l but needed to account for 
the observed characteristics of NGC 2264-C (see end of § 1.2), 
cannot find its origin in protostellar feedback. Although little is 
known about magnetic fields in the region, the most plausible 
remaining hypothesis is that the missing support is magnetic in 
character In any case, gravitational contraction seems to largely 
dominate the physics of the NGC 2264-C clump and may even 
be the main source of the high level of turbulence observed in 
the protocluster (see also Peretto et al. 2006| l. 



5.4. Caveats 

We recall that there is a factor of ~ 8 uncertainty on the out- 
flow momentum flux values given in Table |4] (see § 4.2). This 
uncertainty is mainly due to the uncertain excitation temperature 
of the '^CO(2-l) transition and the unknown inclination angles 
of the outflows. Here, we explore other effects which may have 
led us to underestimate the energy injected by outflows in the 
NGC 2264-C protocluster 

First, our '-CO(2-l) mapping survey may have missed some 
outflows, due to its limited sensitivity. We had a sensitivity limit 
of ~ 0.1 K in our ^^CO{2-\) observations of NGC 2264-C. For 
a typical outflow with a characteristic velocity of ±10 km.s"' 
(with respect to the systemic velocity) and a dynamical timescale 
of 10"* years, this corresponds to a momentum flux sensitiv- 
ity of ~ 0.2 X 10^5 Mo.km.s-i.yr-i. 

A rough estimate indicates that 225 additional outflows as pow- 
erful as the (F1,F2) bipolar flow, or 4 outflows as energetic as 
the LI 551 outflow, would be needed in order to balance gravity 
in NGC 2264-C (see §5.2.2). The sensitivity of our ^^CO{2-\) 
survey is such that we could not have missed such strong out- 
flows. Nevertheless, a large number of fainter, undetected out- 
flows driven by low-mass objects may be present in the proto- 
cluster. 



According to |Peretto et al. (2007), the NGC 2264-C protocluster 
contains seven Class 0-like objects with envelope masses rang- 
ing from 4 Mq to 16 Mq. Based on the Menv-Lboi diagram shown 
in Fig. |9] we estimate that these seven protostellar objects are 
the progenitors of stars with main-sequence masses in the range 
2 Mq < m{""'' < 8 Mq (see § 5.2.1). If flie population of pro- 
tostellar objects in the NGC 2264-C protocluster follows the 
Kroupa (2001 ) IMF, then we expect ~ 250 protostellar objects 
with final stellar masses M{'""' in the range 0.01 Mq i m{'""' * 
2 Mq. Assuming these 250 objects lie in the border zone be- 
tween the conceptual Class and Class I stages and have Menv ~ 
10 of them have envelope masses Menv in 



m{'"''', we infer that 



the range 1 M^ 



Me 



2 Mo, ~ 1 10 of them have 0. 1 M^ 
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Menv i 1 Mq, and -130 of them have 0.01 Mq < Menv ~ 
0.1 Mq in NGC 2264-C. 

Following |Bontemps et al. (19961, we adopt the following re- 
lation between outflow momentum flux and protostellar en- 
veloppe mass: log(Fflow) = -4.15 + l.lx log(Menv)- From this 
relation, we estimate that the ~ 250 low-mass protostellar ob- 
jects mentioned earlier can generate a total momentum flux of 
~3.3xlO"^ Mo.km.s^'.yr"^ Therefore, even if we account for 
the effect of weak, undetected outflows driven by low-mass pro- 
tostellar objects, the total outflow momentum flux injected in 



NGC 2264-C is only F„ 



5x10 Mg.km.s '.yr ', which is 



still a factor of ~ 8 lower than the force needed to support the 
clump against global collapse. 

Second, our '-CO(2-l) map of NGC 2264-C provides only 
a snapshot of outflow activity in the protocluster. Given that the 
characteristic timescale of protostellar evolution is t ~ 10^ years, 
we can estimate the total momentum P,nax injected by protostel- 
lar outflows over their entire lifetime. Bontemps et al. (1996) 
showed that outflow momentum flux is roughly proportional 
to protostellar envelope mass, i.e., Ffiow decreases as Menv de- 
creases and protostars evolve from Class to Class 1 objects. 
Furthermore, the linear correlation Fflo„ oc Menv is suggestive of 
a simple overall accretion/ejection history according to which 
Menv, Fflow, and the mass accretion rate. Mace, all decline expo- 
nentially with time (see Bontemps et al. 1996). We thus com- 
puted the total momentum as: 



-f 

Jo 



F X p"'^'" Ht - F 

A max c UL 1 ry 



X T. 



(14) 



Using the value of the total momentum flux which includes the 
contribution from the low-mass protostellar population, we find 
Pmax - 500 Mo.km.s~*. If we consider a total number N„,ax - 
256 of protostellar objects in NGC 2264-C (including the low- 
mass protostellar population) and adopt a mean stellar mass 
< m* > = 0.2 Mq for the |Kroupa| (2001 ) IMF, then we find 
that the mean momentum per unit stellar mass is: P* = Pmax / 
(NmaxX < >) ~ 10 km.s"'. We conclude that, on average, 
the NGC 2264-C protocluster exhibits a mean outflow momen- 
tum per stellar mass which is only ~ 1/5 of the value Nakamura 
|& Li| ( |2007| l used in their numerical simulations. 

It should be also pointed out that the star formation rate 
within the NGC 2264-C clump may increase in the near fu- 
ture. Using a global mass inflow rate Minf ~ 3x10"^ Mo.yr"' 
for the collapsing clump (Peretto et al. 2006 1, we predict a star 
formation rate SFR = e x Minf ~ 1.5x10 Mg.yr where 
e = 50% is the typical star formation efliciency for individ- 
ual cores ( |Matzner & McKeel|2000| ). Such a star formation rate 
would lead to a momentum injection rate of the corresponding 
outflows SFR X P* ~ 15x10"^ Mo.km.s"^yr"*. This value re- 
mains lower than, but of the same order as, the force needed to 
support the clump against global collapse. This suggests that if 
NGC 2264-C undergoes a burst of star formation in the near fu- 
ture as a result of its global collapse, then feedback due to proto- 
stellar outflows may eventually become sufficient to slow down 
and halt the collapse of the clump. 



1 . We detected a network of eleven CO(2-l ) outflow lobes with 
characteristic projected velocities between ~ 10 km/s and 
~ 30 km/s, lengths between ~ 0.2 pc and ~ 0.8 pc, mo- 
menta between ~ 0.1 Mo.km.s"' and ~ 3 Mo.km.s"', and 
momentum fluxes between ~ 0.5 x 10"^ Mo.km.s^'.yr"' 
~ 50x 10-^ Mo.km.s-i.yr-i. 

2. One pair of high- velocity CO(2-l) lobes corresponds to a 
highly-collimated bipolar flow, whose most likely driving 
source is the millimeter continuum object C-MM13 located 
very close to the center of mass of the NGC 2264-C proto- 
cluster Two other CO(2-l) outflow lobes appear to be driven 
by the millimeter continuum objects C-MMl and C-MM2. 
The driving sources of the seven other outflow lobes have 
not been identified. 

3. The three outflows with identified driving sources have prop- 
erties that are typical of the flows observed toward low- 
to intermediate-mass Class protostars. Based on their po- 
sitions in the Menv - Lboi evolutionary diagram for proto- 
stars, we confirm that the sources C-MMl, C-MM2, C-MM3 
and CMM13 are bona-fide Class objects, which will form 
intermediate-mass stars with m{"'"' ~ 2 - 8 Mq. 

4. The total mechanical power associated with the network of 
eleven ouflows, Ltot = 0.7 + 0.5 L©, is large enough to con- 
tribute a significant, if not dominant, fraction of the turbu- 
lence observed in the NGC 2264-C cluster-forming clump. 

5. However, even if we account for weak, undetected outflows 
driven by low-mass protostellar objects and below the sen- 
sitivity limit of our CO(2-l) survey, the total momentum 
flux generated by the entire outflow population in the proto- 
cluster, Fmax ~ 5x10""' Mo.km.s^^yr^^ falls short by nearly 
an order of magnitude to provide sufficient support against 
global collapse of the clump at the current epoch. Therefore, 
the additional support against gravity missing in the hy- 
drodynamic simulations of NGC 2264-C by Peretto et al 



(2007), cannot arise from protostellar feedback, and most 
likely originates in magnetic fields. 
6. We conclude that gravitational contraction largely dominates 
the physics of the NGC 2264-C clump at the present time and 
may even be the main source of turbulence in the protoclus- 
ter Given the early stage of evolution of the protocluster and 
its global collapse, it is nevertheless plausible that the star 
formation rate will increase in the future, up to a point where 
the feedback due to protostellar outflows becomes sufficient 
to halt the large-scale, global contraction of the protocluster. 
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6. Summary and conclusions 

We carried out a high-resolution '^CO(2-l) survey of 
NGC 2264-C with the IRAM 30 m telescope in an effort to take 
a quantitative census of protostellar outflows in this collapsing 
cluster-forming clump. Additional mapping observations were 
also taken in '^CO(2-l) and C'^0(2-l). Our main results and 
conclusions can be summarized as follows: 
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Table 3. Estimates of the masses of outflow-entrained material at low velocities 



OUTFLOW 


Additional velocity 


Additional mass 


Index of 


Additional 


LOBE 


interval [V2 ; V3] 


^add in [V2 ; V3] 


power-law fit, a 


mass Mio„ m [V3 ; Vsys] 




(km.s ) 


(10 Mo) 




(10 Mo) 




(1) 


(2) 


(3) 


(4) 


Fl (Red) 


[13 ; 10] 


7 ± 2 


-2.7 ± 0.6 


25 


F2 (Blue) 


[1.1 ;4.3] 


1.9 ± 0.5 


-2.5 ± 0.3 


32 


F3 (Blue) 


[1.4 ; 4.9] 


6 ± 3 


-2.3 ± 0.3 


24 


r4 (Ked) 


r 1 "2 . 1 n on 

[13 ; lU.zj 


Z.i ± i.j 






F5 (Blue) 


[0.3 ; 5] 


6 + 2.8 


-3.1 ± 0.9 


31 


F6 (Red) 


[13.1 ; 9.9] 


5.7 ± 3.1 


-3.7 ± 0.8 


26 


F7 (Red) 


[14.7 ; 10] 


15 ±7 


-2.8 ± 0.7 


90 


F8 (Blue) 


[1.7 ; 4.4] 


0.37 ± 0.15 


-2.4 ± 0.4 


7 


F9 (Red) 


[13.7; 9.6] 


1.2 ± 0.7 


-2.4 ± 0.5 


6 


FIO (Red) 


[13.7 ; 9.7] 


1 ± 0.55 


-2.5 ± 0.3 


7 


Fll (Red) 


[15.2; 9.5] 


30± 18 


-2.0 ± 0.9 
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Velocity range [V2 ; V3] in which r;(v)|„/ > 3cr„/ and r;(v)|„„,/to,^. > r;(v)|re/ - So-^/. 

Mass obtained by integration of [T'^{v)\„utfiow - T'^{v)\„f] over [V2 ; V3]. 

Power-law index of the best fit on line wings observed in [Vi ; V3]. 

Estimate of the mass which remains hidden in the main body of the line [V3 ; N syA- 



